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Abstract 

We propose an idea for probing spin states of two coupled quantum dots 
(CQD), by the dc Josephson current flowing through them. This theory 
requires weak coupling between CQD and electrodes, but allows arbitrary 
inter-dot tunnel coupling, intra- and inter- dot Coulomb interactions. We 
find that the Coulomb blockade peaks exhibit a non-monotonous dependence 
on the Zeeman splitting of CQD, which can be understood in terms of the 
Andreev bound states. More importantly, the supercurrent in the Coulomb 
blockade valleys may provide the information of the spin states of CQD: for 
CQD with total electron number N=l,3 (odd), the supercurrent will reverse 
its sign if CQD becomes a magnetic molecule; for CQD with N=2 (even), 
the supercurrent will decrease sharply around the transition between the spin 
singlet and triplet ground states of CQD. 
PACS numbers: 74.50.+r, 73.40.Gk, 73.20.Dx, 72.15.Nj. 
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Coupled quantum dots (CQD) fabricated in semiconductors, or so called "artificial 
molecule", is an ideal quantum system for the study of coherent transport. Each quan- 
tum dot (QD) contains discrete energy levels (quantized due to spacial confinement), which 
can be tuned independently by the gate voltages. The coupling between the dots is also 
adjustable, resulting in either "ionic molecule" for weak coupling or "covalent molecule" for 
strong coupling. The electron-electron interactions play an essential role in the transport 
properties of CQD, due to that the electrons are added to the molecule one by one, and the 
conductance exhibits Coulomb blockade oscillations. Such appealing system has attracted 
a lot of experimental and theoretical attentions, such as the observations of Coulomb block- 
ade in double or triple quantum dots the measurement of microwave spectroscopy 
of a quantum dot molecule and the comparison with the time dependent theory f^-j^; the 
analysis of the additional spectra of CQD in the magnetic field by using single electron 
capacitance spectroscopy P); the investigations of electronic correlations and Kondo effect 
in CQD by two impurity Anderson model P,p!0[|, and much more. 



Nevertheless, the spin -dependent transport through CQD are less addressed |]TT],|T2[. For 
a single QD, it was proposed recently that QD in the Coulomb blockade regime may act 
as an efficient spin filter in the presence of magnetic field [0. For coupled double QDs, 
the spin related physics is even more rich. Consider CQD filled with two excess electrons, 
each dot occupied by one of them because of the large on-site Coulomb repulsion. CQD in 
this case may have four spin configurations: |1 t 2 t), |1 | 2 |), |1 t 2 |), |1 | 2 t), as shown 
in the lower part of Fig.l. Due to the inter-dot tunnehng processes, |1 t 2 |) and |1 | 2 t) 
form a spin singlet state with lower energy. In the presence of magnetic field, however, the 
Zeeman energy may overwhelm the singlet binding energy, and 11121) or |1|2|) becomes 
energetically favorable. Thus, CQD may have either a magnetic or a non-magnetic ground 
state, depending on the external magnetic field. The remaining problem is how to detect 
these states of CQD. We propose that by attaching CQD to two superconductors (S), spin 
states of CQD can be probed by measuring the dc Josephson current. The idea is illustrated 
in the upper part of Fig.l, and hereafter the system is referred to as S-CQD-S. It is well 
known that the dc Josephson current is sensitive to the spin polarization of the weak link 
area of the junction. For the Josephson junction coupled by an Anderson impurity (or S-QD- 
S heterostructure), theories show that the critical current will reverse its sign if the impurity 
(or QD) is singly occupied and therefore becomes a magnetic dot |H|-p^. Recent experiment 
also reported the observation of the vr-junction transition, in a Josephson junction consisting 
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of superconducting banks and a weakly ferromagnetic interlayer . 



Motivated by the above facts and ideas, we will investigate the supercurrent flowing 
through S-CQD-S system. Different from some of the previous works, e.g. [|T^, we take 



f/ ~ A rather than U ^ A, so that the Andreev reflection (AR) process [jl9|,^ and the 
Coulomb blockade (CB) effect may "combine" with each other {U is the constant of the 
intra-dot Coulomb interaction, A is the gap of the superconducting electrodes, and the 
choice of [/ ~ A is physically reasonable). In the calculation, we take into account not only 
intra-dot but also inter-dot Coulomb interactions, both of which are found to be important 
in the fitting of the experimental data 0]. 

The suggested S-CQD-S system can be described by the following Hamiltonian, 

H = Hl + Hdd + Hu + Ht , (1) 

in which 

Hdd = ^ Ei^ni^ + Uini^nii + ^2a^2a + ^^2^21^21 + (2) 

a 

is for the CQD, modelled by two impurities coupled by inter-dot interaction and inter- 
dot tunneling {rii^ = cj^Ci^ is the particle number operator); Hi and are standard 
BCS Hamiltonians for the s-wave superconducting electrodes, with superconducting phases 
4'L = 0/2 and (pR = —(p/2, respectively |2ll]; and Ht is tunneling Hamiltonian , connecting 



the three parts together. 

We begin with the study of the eigen states of isolated CQD. Notice that Hdd can 
be exactly diagonalized in the particle number representation, i.e., in the set of 16 bases 
|0) , |1 I) , ■ ■ ■ 111112 12 1) [|^. Due to the conservation of particle number and the 
conservation of spin in Hdd, the 16x16 space can be divided into several sub-spaces, 

16 = (1)n=o + (2 + 2)n=i + (1 + 1 + 4)n=2 + (2 + 2)n=3 + (1)n=4 • (3) 
We have special interest in the N=2 sub-space, in which Hdd is expressed as 
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with El = En + i?2T + f^i2, E2 = Ei\ +-E21 +f^i2, -£^3 = -^'it +£'21 +U12, E4 = E2t + -£"11 +U- 



'12, 



E. 



Ei^ + Ell + Eq = E21 + E21 + U2- For the case of identical dots p3 



i.e., Ui = U2 = U, U12 = V, E: 



IT 



E. 



2t 



Eq — h, Ell ~ = Eq + h, the eigen 
solution of Hdd has a simple form: {|1 t 2 t), 2Eo - 2h + V}, {|1 i 2 i), 2Eo + 2h + V}, 
2Eo + V}, 2Eo + U}, {|a_), 2Eo + V - E,}, 2Eo + U + E,}, in which 

= -1. (|l t 2 i) T |1 i 2 T)), |/5±) = ^ (|1 T 1 i) T |2 T 2 i)), |a_) and are the two 
new states perturbed from and 11121), |li2|) and are corresponding to 

the S=l triplet states with m=l,-l,0; while |q;_) is the S=0 singlet state, with a binding 
energy 



E, 



^/{U - V f + (4t)2 -{U-V) 



(5) 



With the 16 eigenstates of Hdd, the occupation number of CQD can be evaluated by using 
the relation (o) = Tr{po), in which p = -^e~^^ is the density matrix operator. The left 
and right insets of Fig. 2 shows the occupation number per spin {ria) = (^i<j) + (^20-) vs the 
resonant level Eo for different Zeeman splitting h. (In practice, Eq can be tuned by the gate 
voltage, and h induced by applying an in-plane magnetic field.) The left inset is for the case 
oi h < he, while the right is for h > h^, where 2hc = Eh depicting the competition between 
the Zeeman energy and the singlet binding energy. These curves show that CQD with total 
electron number N=l, 3 is easily magnetized in a magnetic field, while CQD with N=2 
favors a spin singlet state and only transfers to a magnetic state upon a critical magnetic 
field. 

Next, we turn on the weak coupling between S electrodes and CQD. "Weak coupling" 
means that the supercurrent flowing through CQD only serves as a probe to provide the 
information of CQD, without disturbing the quantum states there. To include the tunneling 
between the two dots and the physics of AR, a 4x4 representation is introduced. 
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(6) 



Notice that the retarded Green function of the isolated CQD (g'') can be constructed exactly 
by the Lehmann spectral representation, 

1 



Z ^ UJ - {En- Em) + lO+ 



{m\A\n) {n\B\m) 



(7) 



in which n or m runs over the 16 eigenstates of H £>£)■, and A oi B denotes Qo- or cj^. 
We assume that the full Green function of S-CQD-S (G*") can be derived by the following 



approximation 24 



(8) 



in which S*" is the self-energy caused by the coupling between S electrodes and CQD. S*" is 
obtained as 



'^S^ ^ 
S 



(9) 



in which 
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(Im^/x > 0) 



(10) 



(11) 



and r^/zj is the coupling strength between left / right electrode and CQD. Then the dc 
Josephson current flowing through S-CQD-S is obtained by the Green function technique as 



/ = y sin0 j ^f{uj)j{uj) 



in which 



j{uj) = Im J(u;) = Im 



(^ + iO+)2_A2 det[(g 
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(12) 



(13) 



and f{uj) = l/{e^^ + 1) is the Fermi distribution function. det[(g'^)~^— H''] in the denom- 
inator has several real roots (with infinitesimal imaginary part -iO+) within the range of 
\uj\ < A, corresponding to the Andreev bound states (ABS) in the S-CQD-S system. To 
avoid the divergence around these singularities, we adapt the integral path to a V-shaped 
contour shown in the middle inset of Fig.2 [ES|1, result 

/ = ^sin01m / ^/MJM . (14) 

n Jy ZTT 

Fig.2 shows the critical current Ic = I {(f) = f ) vs i?o with different h. The curve of h = 0, 
as expected, has four CB peaks, among which are the valleys of the total electron number 
from N=0 to N=4. ( The curve is symmetric to Eq = —(U/2 + 2V) by virtue of the electron 
and hole symmetry.) With the increase of h, the 1st and 4th peak are gradually suppressed 
|26|, while the 2nd and 3rd peak exhibit a non- monotonous dependence on h: decrease first 
and reverse its sign to a negative peak, then increase again and diminish at sufficient large 
h. To understand this anomalous h dependence, we plot the curves of Ic vs h in Fig.3a, for 
a non-interacting S-CQD-S system by setting U = V = 0. (In this case, the approximate 
Eq.(8) becomes an exact one, and g*" can be evaluated explicitly.) For h = 0, there are two 
peaks in the curve separated by 2t; for h > 0, Ic a-t Eq = also exhibits a non-monotonous 
h dependence (see the up-right inset). As we know, the supercurrent fiowing through CQD 
is conducted mainly by the ABS, which have the property that two adjacent states carry 
the supercurrent with opposite signs. The spectrums of j(u;) for h = 0,t,2t are shown in 
the insets of Fig. 3a. In each spectrum, there are four ABS denoted by 1"^ and 2^ within 
the superconducting gap, and continuous spectrum outside the gap. For /i = 0, 1^ and 
2^ are distributed symmetrically to the Fermi surface fi = 0. With the increase of h , these 
states move down toward u = —A. The case of h = t corresponds to three of ABS 1^ 
and 2~ are below the Fermi surface, while h = 2t corresponds to all of them below the 
Fermi surface. By taking account of the Fermi distribution and the continuous spectrum 
contribution from c^, the anomalous h dependence is readily understood. Considering the 
Coulomb interactions will induce more ABS in ji^uj), but the h dependence of the peaks in 
the Ic vs Eq curve can be explained in the similar way. 

We are more interested in the supercurrent fiowing in the CB valleys, since the spin states 
of CQD are well defined there. The curves of Ic vs h and corresponding m = (n^) — (rii) vs h 
are shown in Fig. 3b, with Eq chosen in the CB valleys of N=l, 2, 3. The two curves marked 
with A are typical for the supercurrent fiowing in the valley with odd number, where CQD 
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has a net spin. Due to the strong Coulomb interaction, small magnetic field {h ~ ksT) will 
lead to the transition of CQD from a non-magnetic molecule to a magnetic one; meanwhile, 
Ic reverse its sign and experiences a "tt-J unction" transition. The two curves marked with B, 
in contrast, are typical for the supercurrent fiowing in the valley with even number, where 
CQD may choose a non-magnetic singlet or a magnetic triplet as its ground state. The 
transition between them occurs at the critical value of Zeeman splitting 2hc = Eb, with 
Ic > for h < he and Ic — > 0"*" for h > he- Out of our expectation, for h ^ he, le —>■ 0"*" 
instead of /c ^ 0~. Roughly, this might be interpreted as the addition of two 7r-junction 
transitions, namely, the supercurrent reverse its sign after fiowing through the first magnetic 
dot, but returns back after fiowing through the second one. 

To sum up, we have investigated the dc Josephson current fiowing through S-CQD-S 
hybrid system. Our calculation is based on the exact diagonalization of coupled two impurity 
model and the construction of Green function in the spectral representation. To avoid the 
difficulty of solving the exact ABS, we choose a V-shaped integral contour in the complex 
uj plane. In the numerical study, wc find that the dc Josephson current flowing through S- 
CQD-S can provide rich information of the spin polarization of CQD. The CB peaks exhibit 
a non-monotonous dependence on the Zeeman splitting of CQD, the current in the CB 
valleys can be used to probe the spin states of CQD. We believe that the proposed S-CQD-S 
system is within the scope of the update nano-technology in S/2DEG heterostructure, and 
we are looking forward to seeing the relevant experiments. 

This project was supported by NSFC under Grant No. 10074001. 
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FIGURE CAPTIONS 



Fig. 1 Upper part: Schematic diagram of the proposed S-CQD-S system. Lower part: The 
four spin configurations of CQD with the total electron number N=2. (a) and (b) are 
for |1 t 2 t) and |1 J, 2 J,), respectively; (c) illustrates that |1 f 2 J,) and |1 J, 2 f) are 
coupled via high energy virtual states and forms a spin singlet. 

Fig. 2 The critical current /c = /(|) vs the resonant level Eq with different Zeeman sphtting 
h, in units oi e — h — A — 1. The parameters of CQD are: U — 0.6, V — 0.2, 
t — 0.1, due to which /ic ~ 0.04. The temperature and coupling strength are chosen as 
ksT — 0.01 » r = 0.001, as required by the weak coupling limit. The solid, dashed, 
dotted, and dash-dotted curves correspond to h — 0.00, 0.02, 0.06, 0.10, respectively. 
The left and right insets show the occupation number per spin (n^r) vs Eq for h = 0.02 
and h = 0.06, respectively. The middle inset schematically shows the singularities of 
J(a;) and /(a;) and the V-shaped integral contour. 

Fig. 3 (a) Ic vs Eq curves with different h, for the case of non-interacting S-CQD-S in 
which U = V = 0. The solid, dashed, dotted, and dash-dotted curves correspond 
to h = 0.00, 0.05, 0.10, 0.20, respectively. Other parameters arc the same as those of 
Fig. 2. The up-right inset shows vs h at Eq = 0, which exhibits a non-monotonous 
dependence. The other three insets show the spectra of j(a;) at Eq = 0, with different 
h marked in the plots, (b) Ic vs h (solid) and corresponding m = (n-f) — (nj^) vs h 
(dotted), in the Coulomb blockade valleys of N=l, 2, 3. The two curves marked with A 
are for Eq — —0.95 (also Eq — —0.05), typical for CQD with odd number of electrons; 
the other two marked with B are for Eq — —0.50, typical for CQD with even number 
of electrons. 
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